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This work considers the effects of the Hurst exponent (H) on the local electric field distribution 
and the slope of the Fowler-Nordheim (FN) plot when considering the cold field electron emission 
properties of rough Large-Area Conducting Field Emitter Surfaces (LACFESs). A LACFES is 
represented by a self-affine Weierstrass-Mandelbrot function in a given spatial direction. For 0.1 ^ 

H < 0.5, the local electric field distribution exhibits two clear exponential regimes. Moreover, a 
scaling between the macroscopic current density ( Jm ) and the characteristic kernel current density 
(Jkc ), Jm ~ [Jkc]^ H , with an H-dependent exponent /3h > 1, has been found. This feature, 
which is less pronounced (but not absent) in the range where more smooth surfaces have been 
found (0.5 ^ H ^ 0.9), is a consequence of the dependency between the area efficiency of emission 
of a LACFES and the macroscopic electric field, which is often neglected in the interpretation of 
cold field electron emission experiments. Considering the recent developments in orthodox field 
emission theory, we show that the exponent /3h must be considered when calculating the slope 
characterization parameter (SCP) and thus provides a relevant method of more precisely extracting 
the characteristic field enhancement factor from the slope of the FN plot. 
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Introduction 

The theory of cold field electron emission (CFE) is relevant for understanding the important technological as¬ 
pects in the development of large-area electron emitter devices with apex radii of no less than a dozen nanometers 
[I] . This theory can be conveniently formulated using a free-electron model for the emitter and a Schottky-Nordheim 
barrier model for the tunneling barrier m- Large-area field emitters provide an effective scenario for complying with 
the requirement that electron field emission in metals occurs when a large local electric field (~ a few V/nrn) induces 
electron tunneling through a potential barrier out of the corresponding surface into vacuum. Therefore, practical 
research on CFE involves the development of techniques to reduce the barrier through which the electrons must 
tunnel [3i. For example, an interesting phenomenon of electron emission at low macroscopic electric fields has been 
reported for many materials, such as amorphous diamond-like films [7]. Diamond cold cathodes formed using chemical 
vapor deposition (CVD) for use in field emission displays have attracted considerable attention because of their low 
work function mm]. Moreover, the local field distributions over emitting surfaces have been explored, for example, 
using carbon-nanofiber (CN) films constructed of highly uniform nanofiber arrays; in that case, experimental evidence 
of Gaussian behavior of the distribution of the spatially resolved field enhancement factor was obtained via scanning 
anode field emission microscopy (SAFEM) [TO] . In addition, experimental results obtained using the same technique 
(measuring the constant macroscopic field emission current map) have revealed exponential decay in thin films of 
non-oriented carbon nanotubes (CNTs), which exhibit an irregular morphology on the corresponding surfaces |llj . 
Indeed, it is becoming increasingly clear that modified Fowler-Nordheim-type (FN-type) equations must be developed 
to address CFE from rough thin-film emitters (for example, films of CNTs and related materials) |12) . In such a case, 
upon the application of a voltage to the support tips, a high local electric field develops at the apex, inducing a field 
emission tunneling current m • In fact, no definitive treatment has yet emerged for CFE from conducting materials 
that accounts for roughness at small scales. In these circumstances, improving the Fowler-Nordheim-type equation 
to yield a satisfactory fit to the experimental results is an important requirement for the development of field emitter 
devices with more precise specifications. 

Only a relatively small number of studies have considered the role of the irregular morphology of the surface of a thin 
film on the local electric field distribution HMS3 and the corresponding field electron emission properties of Large- 
Area Conducting Field Emitter Surfaces (LACFESs) in a genuine three-dimensional problem [16, .17] ■ However, 
considering the morphologies of the conducting surfaces that are commonly experimentally investigated, a large 
number of problems exist in which such a consideration is relevant usms]. Irregular morphological features that 
remain on the surface after the manufacturing process may act as field emitting tips in the presence of an external 
electric field. Indeed, a more detailed investigation of the surface preparation prior to the application of a high electric 
field leads to a considerable decrease in the breakdown rate [20 . 

Motived by the aforementioned studies, in this work, we present the first consideration of the role of the Hurst expo¬ 
nent H of an irregular LACFES in determining the related local electric field distribution and field electron emission 
properties using orthodox field emission theory. A LACFES is represented by a self-affine Weierstrass-Mandelbrot 
(WM) function in a given spatial direction such that all surfaces are assumed (i) to exhibit no roughness along any 
other spatial direction and (ii) to have the same global roughness (which is a measure of the fluctuations at large 
scales). The latter condition allows us to systematically evaluate the role of H (and the short-wavelength fluctua¬ 
tions in the morphology of the LACFES) in determining the non-linearities that may appear in the corresponding 
Fowler-Nordheim (FN) plots. 

The results indicate that for any H, the local electric field distribution over a LACFES exhibits two exponential 
regimes, implying a non-linear (power-law) relation between the macroscopic current density and the characteristic 
current density, if orthodox field electron emission is assumed. We show that this scaling is a consequence of the 
dependence between the area efficiency of emission of the LACFES and the macroscopic electric field, which is more 
pronounced for 0.1 ^ H < 0.5. This result allows for the introduction of a new slope correction parameter in an 
FN-type equation to allow for the more precise extraction of the characteristic field enhancement factor (FEF) from 
the slope of the FN plot. 

Morphology of a LACFES 

Central to our approach is an expression that represents the height profile of the emitter along a line on the 
emitter surface (defined here as the x coordinate). The other coordinate parallel to the emitter surface is denoted by 
y , and the coordinate normal to the surface is denoted by z. To model this profile, we use the Weierstrass-Mandelbrot 
(WM) function [2l[ : 
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FIG. 1. Scale of the local roughness for the WM function with H £ [0.1, 0.7]. The very good agreement between the slope ( ai ) 
of log 10 [u;(r)] as a function of log 10 [r] and the parameter H is remarkable. From top to bottom, the dashed lines correspond to 
slopes of ai = 0.1,0.2,0.3,0.4, 0.5 and 0.7. In the inset, a Large-Area Conducting Field Emitter Surface (LACFES) represented 
by a self-affine WM function (see Eq. 0 ) with the parameter H = 0.5 and £ = 0.5e is shown, where “e” is Neper’s number. It 
is evident that the conducting surface lies on a two-dimensional substrate with no roughness along the y direction. 


N 

}wm(x) = cY^- Hn sin(K 0 t; n x + (1) 

n—0 

where c is a constant related to the global roughness amplitude, 0 < H < 1 is the Hurst exponent parameter (also 
called the “local roughness exponent” (cq) in experimental thin-film science), Kg is the fundamental wavenumber, 
and (f) n is an arbitrary phase (in this work, we consider the value of 0 < cf> n < 2 tt to be randomly chosen). £ is a 
parameter (£ > 1) which is a measure of the distance between the frequencies. This function can be used to model 
the emitter height profile as a superposition of a set of sine waves, with the wave amplitude decreasing as £ increases. 

Larger values of the Hurst exponent imply a smoother surface because the ratio of the logarithm of the local 
roughness (Eq. Q below) to the logarithm of the scale length is larger (here, the scale length refers to the scale at 
which the statistical properties, such as height fluctuations, of an irregular surface are analyzed; i.e., the box size “r” 
as a measure of the local roughness is defined in Eq. 0. For 0.5 < H < 1, the behavior is called persistent, whereas 
for 0 < H < 0.5, the behavior is called is antipersistent. Antipersistence indicates that the heights of various points 
on the surface are negatively correlated and that the correlations rapidly decay to zero. This scenario corresponds to 
a small ratio of the local roughness with respect to the scale. Finally, at H = 0.5, the heights of the LACFES are 
statistically random. If we use the one-dimensional correlation function Q(x), defined as 


fl(x) 


-fwM{~x)fwM{x) 
[, fwM{x )] 2 


( 2 ) 


where the overbars represent spatial averages, it is trivial to show that fl(x) does not depend on x. In fact, using 
Eqs. (0 and 0,0= — l). Thus, the signal of 0 defines the type of correlations that appear in the WM 

function, in such that O > 0, O < 0 and O = 0, for H > 0.5, H < 0.5 and H = 0.5, respectively. The WM function is 
a combination of periodic functions that exhibits two major features: (i) it can be anisotropic, and (ii) it is self-affine 
with complexity at small scales. Thus, the consideration of such surfaces is motivated by the fact that real irregular 
surfaces are generally neither purely random nor purely periodic. 

Now we use the notation h = /wm(x ) to represents the height of the profile in the 2 direction with respect to the 
plane h(x,y) = 0. The global roughness parameter ( W ) provides a measure of the height variation in the emitter 
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FIG. 2. Profile of a LACFES with H=0.5 (shown in Fig. |T]) and equipotential lines calculated numerically from the solution 
of the Laplace equation with appropriate Dirichlet conditions (4> s = 0 V and <F 4 = 500 V - see the text for more details). In 
the inset, a portion of the LACFES and the equipotential line defined by = 1 V, representing the approximation of the 
field emitter surface, are highlighted. The horizontal dashed (red) line indicates the average height of the roughness profile. 
The far-away conducting anode (not shown) is located at h ss 16/xm from the average height of the LACFES (see the text for 
more details). 


profile, taken across the emitter as a whole. This parameter is defined as the root mean square of the front-surface 
height fluctuations and is given by 


(h-ti) 2 


(3) 


where the averaging is performed over a length L x that represents the lateral size of the emitter in the x direction. 
To help determine whether the WM function provides a realistic and useful description of self-affine one-dimensional 
profiles acquired from real emitter surfaces, we can analyze how the local roughness u>(r) scales with the length r, which 
is a useful measurement for experimental purposes. This is also defined as a root-mean-square height fluctuation, but 
the spatial average is limited to a scale of size r, with r < L x . More explicitly, the local roughness is given by 


u(r) = ([(h-h) 2 ]V 2 ), (4) 

where the angular brackets represent the configurational average obtained as the box (of length r) scans the entire 
irregular surface. 

The function fwM{x) has no characteristic length (beyond the size of the system itself), in the sense that the level 
of detail of this function is self-similar under an affine scaling (in which the “x” axis is stretched by a factor of £ and 
f wm(x ) is modified by £ H , i.e., /wm(^x) = £ H /wm(x))- Therefore, its statistical properties are identical at different 
scale lengths. As indicated by the analysis presented in Fig. [T] for small scale lengths, the local roughness scales as 
uj(r) ~ r“ ! , which is in contrast with large scale lengths, for which the local roughness coincides with the global W. 
The very good agreement between the slope cq (from log 10 [w(r)] as a function of log 10 [r]) and the parameter H can 
be clearly observed. In the inset of Fig. [I] the surface that mimics the LACFES in the case of H = 0.5 is shown; this 
surface was computed using N = 200 in the sum defined in Eq. |T]). We stress that x is considered to be an integer 
number that represents a distance measured in terms of some basic unit distance u, as are y and 2 . 

In fact, the typical Fowler-Nordheim-type theory of CFE was not developed at the atomic scale and thus does 
not apply to emitters that are sharp on that scale. To overcome this difficulty, we developed an algorithm to round 
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off the sharpest projections and replace the real WM surface with an equipotential surface, (very similar to the 
previous one), that is calculated from the numerical solution of the Laplace equation using a finite-difference scheme 
(Liebmann method - see Supplementary Information I). For practical applications, this approximation corresponds, 
for instance, to the condition in which good Spindt arrays [22 ] are produced by slowly increasing the voltage such 
that the sharpest tips “burn off.” Our numerical solution has been shown to yield an FEF for ideal protuberances 
that is in very good agreement with the results of previous analytical, finite-element and multigrid methods [23H251 . 

In Fig. [2| several equipotential lines are shown that consider electric potentials for the LACFES (for the case of 
H = 0.5) and the far-away conducting anode (located at h ss 16/rrn from the average height of the LACFES) of 
d 15 = 0 V and = 500 V, respectively, corresponding to Dirichlet conditions. The equipotential surface, on which 
the local electric field distribution will be calculated, is shown in the inset of Fig. [2] In this work, we use u = 5 nm, 
which results in a global roughness of W ss 75 nm for all LACFESs. Moreover, this methodology ensures that for 
small scales (though larger than the atomic scale, e.g., r = 20 nm), the local roughness is negligible, and for large 
scales (i.e., r > 500 nm, although still considerably smaller than the lateral size of the system L x ), the local rough¬ 
ness scales with the same exponent H. The same procedure is adopted for any H that allows for the use of CFE theory. 

Orthodox CFE Theory - Local emission current density of a particular lateral loca¬ 
tion on a LACFES 

The recent developments in CFE theory can be viewed as a relevant approach for understanding the field emis¬ 
sion properties of a rough LACFES. The local emission current density (LECD), A, in the CFE regime can be 
written as a function of a convenient set of parameters, which includes the sum of all traveling states incident on the 
inside of the emitter surface considering all occupation and transmission probabilities. The general result takes the 
following form [25] : 


Jl = ZD , (5) 

where Z is an effective incident current density and D is a transmission probability. In the case of high electric fields 
and low temperatures, the value of D at a specific forward energy level (often that for a Fermi-level electron moving 
normal to the emitter surface) is chosen, and Z is calculated as Jl/D. CFE is an LECD regime in which most 
electrons escape via deep tunneling (i.e., tunneling well below the top of the barrier) from states close to the emitter 
Fermi level, and the Landau and Lifschitz approximation can be applied [26] , Then, using this approximation, Eq. 
([ 5 ]) can be written as follows: 


Jl = Z F P F exp [-Gf], (6 ) 

where “F” is used to indicate parameters related to a barrier of zero-field height (</>). G F is the related barrier strength 
(also called the Gamow exponent or the JWKB (Jeffreys-Wentzel-Kramers-Brillouin) exponent) and is given by 


Gf 


z)dz = g e 


M 1 / 2 {z)dz. 


(7) 


In Eq. ([ 7 ]), ip 2 = 2m e /H 2 , where m e is the electron mass and h is related to Planck’s constant; M(z) specifies the 
form of the barrier, which is the difference between the total electron potential energy, U(z), and the total energy 
component, E z , related to the motion of the electron in the z direction (perpendicular to the surface of the emitter in 
the case of a conducting emitter), M(z) = U(z) — E z . Moreover, the integral in Eq. 0 is performed over the range 
of z for which M(z ) > 0. 

Considering metal emitters with appropriate dimensions, as previously discussed, the LECD J F F) is formally 
given by the following equations: 


Jl{0,F) = \c F EJk(<t>,F ), 


( 8 ) 


J k ((f>,F) = acj) 1 F 2 exp 


—v F b(j) 3 / 2 /F 


(9) 


where </*,(</>, E) is the kernel current density, which can be evaluated exactly for any chosen barrier form, i/p, and 
choice of <t> and F; a (=1.541434 x 10~ 6 A eV V~ 2 ) and b (=6.830890 eV ~ 3 / 2 V nm -1 ) are the first and second 
Fowler-Nordheim constants, respectively [25] , 
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In this work, we consider the Sclrottky-Nordheim (SN) barrier used in the Murphy-Good theory [271 and, more 
recently, in the Forbes approximation (orthodox emission theory) [2] . The SN barrier corresponds to the lowering of 
the top of the barrier relative to that of an exactly triangular barrier such that M(z) = <f> — eFz — e 2 /16ire 0 z (e is 
the positive elementary charge and eo is the electric constant). In this way, it is possible to define the scaled barrier 
parameter, /, which is given by 

/ = F/F r = (4tt e o /e 3 0 2 )F. (10) 

The reference field F R is the field that is required to lower the barrier height by an amount equal to the local work 
function <fi. For any tj>, it is possible to define a parameter p{(f>) as follows: 


?7(0) = b(j)^ 2 /F R » 9.8362(eL/^) 1 / 2 . (11) 

If we compare Eqs. ([6]) and (|8j), the pre-factor Pp used in the first equation is included in A cfe- One problem with 
attempting to obtain good predictions of the CFE current density is that exact values of A cfe (which depends on 
the material) are not well known. The current best guess is that A cfe lies in the range 0.005 < Xcfe < H [28] . 
We compute the local current density that results in A cfe = 1 in Eq. ^ such that Jp(^>,F) = Jk(<j>,F). In any 
case, the correction factor with the largest influence on Jp is the barrier shape correction factor up, which will also 
be considered in this work. Finally, we use an approximation such that each point on the LACFES represents a 
particular lateral location on the emitter surface. 

Results and discussion 

Local electric field distributions 

Before discussing the results for the LECD, we analyze the behavior of the local electric field intensity distribu¬ 
tion over the LACFES, p(F). In Fig. [3] (a), the behavior of p(F) is shown for a LACFES with H = 0.1 and for several 
values of the anode electric potential $" 4 . The vertical black dashed line (gray dashed-dotted line) represents an 
electric field intensity of 2.5 V/nm (4 V/nm - see Ref. [23]), which corresponds to a typical value for field emission in 
pure metals with a local work function of approximately 3.5 eV (4.5 eV). The inset of this figure shows the behavior 
of the mean electric field (F) over the LACFES, which is defined by 

(F) = — FjjAxjAxj , (12) 

A-m 

t-,3 

as a function of <I>‘ 4 for FI £ [0.1,0.9]. Here, Am is the total “substrate footprint” area, F i:j is the local electric 
field intensity at a particular location (*, j) on the LACFES, and Ax,Ax : j = u 2 corresponds to an area of unity, with 
A or,;, Axj <C L x . The results clearly suggest a linear dependence whose slope, as a non-linear function, depends on 
the H of the LACFES. It is interesting to observe the approximate collapse of the (F) behavior for 0.5 ^ FI ^ 0.9. 
This allows us to more easily identify the primary differences between the electrical properties of a rough correlated 
(and random) LACFES (0.5 ^ H ^ 0.9) and a rough anti-correlated LACFES (0.1 ^ IF < 0.5). 

As a measure of emitter sharpness, we compute the distribution of the macroscopic field enhancement factors 
(FEFs), 7 = F t j /F m . that are typically applied in measurements using scanning anode field emission microscopy. A 
high FEF factor leads to a low turn-on voltage and a high emission current, which are desirable for emitter applications. 
In these calculations, the field at the position of the average height of the LACFES (ss 0.37/rm for a LACFES with 
any H - see Fig. [5] for = 0.5), in the absence of protrusions, was considered to be the macroscopic electric field Fm- 
Thus, the flat anode, which is placed at a height of h A (measured with respect to the substrate on which the film 
was grown, i.e., the plane h(x,y) = 0), is located at a distance of (d) = h A — h from the cathode. The macroscopic 
electric field, Fm, is given by 


Fu = W 


(13) 


If <F‘ 4 is taken to be 500 V and (d) ~ 16/im (see Fig. [2]), then the macroscopic electric field is F M ~ 31 V//xm. 
In Fig. [3](b), the results for the distribution of 7, p( 7), are shown for H £ [0.1,0.9]. This measurement can be 
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H 

Sh x 

logio(e) 

Sh x 

l°gio( e ) 

7 c 

0.1 

0.343 

± 

0.001 

0.543 

± 

0.001 

K, 11.86 

0.2 

0.311 

± 

0.001 

0.685 

± 

0.006 

« 10.35 

0.3 

0.315 

± 

0.002 

0.744 

± 

0.006 

« 9.34 

0.4 

0.333 

± 

0.004 

0.93 

± 

0.02 

« 8.95 

0.5 

0.321 

± 

0.005 

1.09 

± 

0.01 

« 8.51 

0.6 

0.317 

± 

0.004 

1.15 

± 

0.01 

« 8.44 

0.7 

0.319 

± 

0.002 

1.24 

± 

0.02 

« 8.12 

0.9 

0.311 

± 

0.005 

1.3 

± 

0.1 

« 7.78 


TABLE I. Results of the extraction of parameters from the p( 7 ) distribution shown in Fig. [ 3 ] (b) for 0.1 ^ H ^ 0.9. For 
0.5 ^ H ^ 0.9, a dominant exponential decay can be observed in the p( 7 ) distribution, characterized by 6 { r (see the text for 
more details). 


regarded as a probability distribution because the number of field emitter locations per unit area on the LACFES 
can be written as dN = p(^)dx, where dN is the number of locations with FEF 7 € [7, 7 + d"f] per unit area. For a 
given value of H, the corresponding distributions p( 7) were found to be the same for any d 1 " 4 , indicating that these 
results are not dependent on $ A . This is the result that is physically expected because in classical electrostatics (in 
the case of a diode), the value of the field enhancement factor depends only on the electrode geometry and not on 
the applied electrostatic potential difference. However, it is a useful test to confirm that this expected physical result 
is reproduced in our numerical simulations. An interesting alternative to this behavior is to observe a collapse of the 
local electric field intensity distributions p(F), which is expected to occur if the variable F is replaced by a scaled 
variable x = F/(F). This phenomenon can be clearly observed in Fig. [3jc), which presents the results for H = 0.1 
and for all <& A values shown in Fig. 0 a )- This collapse was also verified for all H values considered in this work. 

Interestingly, for all H , two exponential decaying regimes are apparent, labeled as (I) and (II), which correspond 
to p T ( 7) ~ exp(— Sjjj) and p n ( 7) ~ exp(—<5^/7), respectively, where SjJ > 1 5 ! H , in general. In particular, for H = 0.1, 
5q A x log 10 (e) = (0.343± 0.001) is identified for interval (I) and 5^ x log 10 (e) = (0.54±0.01) is found for interval (II) 
(the slopes of the dashed-dotted lines in Fig. [30), where “e” is Neper’s number). In fact, AS = 5\j — S T H assumes the 
lowest value for H = 0.1 (see Table I), suggesting that the area in which electron emission occurs is larger in the limit 
of large values of This result, which is a consequence of the first (and second) slower decay of p{ 7), suggests a 
greater probability of finding emitting sites on LACFESs (primarily for H = 0.1) with larger FEFs. We also compute 
the characteristic FEF, 7 c = Fq/Fm , where Fq = max{Fij}. For H = 0.1, it is found that 7 c ~ 11.86. The same 
behavior of the two exponential regimes is observed for all H. 

Table I summarizes the values of the parameters extracted from the ,0(7) distributions presented in Fig. [3] (b) and the 
7c values for all values of H explored in this work. Interestingly, these features capture important experimental results 
for thin-film emitters with irregular surfaces, in which p(py) has been shown to exhibit an exponential dependence on 
the FEF [i30]. However, we note that the FEF values reported in that work were typically between 100 and 300 (i.e., 
approximately one order of magnitude larger than those found in this work). This effect can be understood based on 
the following considerations: (i) the differences between the global roughness of the experimental surfaces used in the 
previous experiment and ours and (ii) the one-dimensional field variation across the emitter surface in our approach 
(see Supplementary Information II). In the cited work [30], the global roughness on a length scale of 20 pm was found 
to be approximately W « 4pm, whereas in our work, the considered LACFESs have a roughness of W ss 75 nm (i.e., 
~ 50 times smaller) on the same length scale for all H (see Fig. [l]) . 

In Supplementary Information II, we present strong evidence to corroborate speculation (ii) by simulating a genuine 
three-dimensional surface using a possible experimental setup modeled as an irregular emitting surface iteratively 
generated by a fractional Brownian motion ( FBM) algorithm with H = 0.1. In this case, again, an exponentially 
decaying regime in p{ 7) is observed. This finding indicates that our LACFES model captures the main features 
of the experimental results, namely, the local electric field distribution and the field emission properties (in the 
orthodox theory) of irregular conducting surfaces. This discussion suggests the importance of correctly measuring 
the H exponents and the global roughnesses of the rough surfaces used to experimentally represent LACFESs (such 
as randomly oriented carbon thin-film emitters) in terms of, for instance, the local roughness scale (see Eq. ©) 
extracted using a probe microscopy technique such as electrostatic force microscopy (EFM). In such a case, a careful 
analysis must be performed to account for the distortions in the image that may occur due to the finite size of the 
EFM tip [3XJ, which may result in overestimation of the Hurst exponent (or the local roughness). 
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FIG. 3. (a) Local electric field intensity distribution, p(F), for a LACFES with H = 0.1. The vertical (black) dashed and (gray) 
dashed-dotted lines represent electric field intensities of 2.5 V/nm and 4 V/nm (see Ref. [29]), which correspond to typical 
values for field emission in pure metals with local work f unct ions of approximately 3.5 eV and 4.5 eV, respectively. The inset 
shows the behavior of the mean electric field (F) (see Eq. ( | 12| ) ) as a function of the anode electric potential <t> A for Ft = 0.1 (red 
circles), H = 0.2 (blue triangles), H = 0.3 (green squares), Ft = 0.5 (black asterisks), H = 0.6 (magenta diamonds), H = 0.7 
(gray pentagons) and F[ = 0.9 (dark yellow right triangles), (b) Local FEF distribution, p(y), for FI £ [0.1, 0.9]. For a LACFES 
with 0.1 ^ 77 < 0.5 (top panel), two pronounced exponentially decaying regimes, (I) and (II), can be clearly observed, and 
the dashed-dotted lines illustrate this behavior for Ft = 0.1. For a LACFES with 0.5 ^ Ft ^ 0.9 (bottom panel), it is evident 
that (I) and (II) correspond to slow and fast exponential decays, respectively, with the latter reflecting the low probability of 
finding high peaks on the corresponding LACFES. (c) Distributions of x = F/(F) for Ft = 0.1 for all <f> A values considered in 
(a). The collapse of the curves is evident, indicating that the results for the p(q) distributions (for all Ft) do not depend on 
the anode electric potential $‘ 4 . 


Area efficiency of emission and Fowler-Nordheim plots 

Previously, the occurrence of a relatively slow decay (interval (I)) was observed in the p(y) distribution for ar¬ 
bitrary Ft, as shown in Table (I). Remarkably, this tendency appears to be characterized by approximately the same 
parameter Sjj for any Ft in this interval. This interesting result suggests that the area in which electron emission 
occurs in a LACFES may be higher compared with the ideal metal surface morphologies that are often modeled 
to explain experimental results. If a metal has no appreciable irregularities on the nanometer scale, this quantity 
corresponds to the Area Efficiency of Emission (AEE), a«, which is typically less than 10 -5 (see Ref. [4]). This 
indicates that for a LACFES with a regular array of single emitters, it is expected that the effective emission area (the 
emitting area at the tips) will be considerably less than the apparent “macroscopic” geometrical area (or “substrate 
footprint” area Am) of the physical emitter that is observed visually. Based on these arguments, if the emission 
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area is, in fact, constant, it is convenient to define a« as follows (here, um = Am, using Forbes’ notation (see Ref. 
0)> because we set A cfe = 1; A m is the macroscopic “pre-exponential correction factor” or, alternatively, “formal 
emission efficiency”): 


_ Jm 

OLM = -, 

JkC 


(14) 


where the characteristic kernel current density , J k c (which is equal to the characteristic local current density, Jc, 
using Forbes’ notation, by the same arguments presented above), and the macroscopic current density, Jm , for a 
LACFES are given by 


JkC = max{J k {4>, F)}, 

where J k ((f>,F) is given by Eq. Q and, because A cfe is approximated as equal to 1, 


(15) 


y Jk(zp)^Ap 

Jm = i/A M = —--• (16) 

■A-m 

Here, “i” is the total emission current, and the sum is taken over the “substrate footprint” area of the emitter, Am- 
Furthermore, the values of J k (z p ) are given by Eq. |9]). These assumptions are equivalent to assuming that the 
empirical CFE (i,^) characteristics should obey 02] 


i = C($ j4 ) k exp (—7V/<!>‘ 4 ), 


(17) 


where C, N and k are constants. If the electron field emission is orthodox, considering the SN barrier, and the 
emission area is constant, then n is expected to be n = 2 — rj(<j))/6, where ??(</) is given by Eq. (11) (see Ref. [32]). 
In Eq. ( [ 17 ] ), we assume that the emission quantities and the measured quantities (*,<1> A ) are identical. 

Note that the linear dependence between Jm and J k c for a rough LACFES is not clear in our case because the 
LACFES exhibits irregularities on small scales (with the apex radii of the emitters greater than dozens of nanometers). 
This question, for the geometries used in this work, has not been previously addressed; therefore, we must address 
it. Studies that consider smooth surfaces generally follow the implicit assumption that om is only weakly field 
dependent, meaning that for practical purposes, one can take it to be nearly constant. However, we show that this 
assumption does not hold, particularly for LACFESs with high local roughness (0.1 < H < 0.5). As we will discuss, 
the consequences of the scale relation between Jm and J k c can provide relevant information concerning the effects of 
the morphologies of rough experimental LACFESs on the FN plot slopes if the electron field emission is orthodox. 

We emphasize that we use an approximation in which each point on the cathode surface represents a particular 
lateral location on the emitter surface, as previously discussed. Moreover, in our calculations, space is filled with a 
discretized cubic lattice with a unit volume of Ax x Ay x A z. Thus, Eq. (16 1 can be written as 


y JL(z p )Ax p Ay p 


Jm 


(N x - l)Ax x (N y - 1)A y’ 


(18) 


where N x and N y represent the numbers of points defining the LACFESs in the x and y directions, respectively. In 
Eq.([T8|, AxiAyi = ... = Ax Nx Ay Ny = AxAy. 

For LACFESs (for all H ), we follow the definition of Eq. (15) by identifying “C” with the apexes of the more 


prominent emission sites. In Eq. (16), if we use the SN barrier of unreduced height </> and the Forbes approximation 
[2], then vp = v(J) ~ 1 — /+ (l/ojfln(f), with /given by Eq. (10). We assume that for the LACFESs considered in 
this work, the local work function is approximately constant over the surface, and we adopt (f> = 3.5 eV. Therefore, 
the parameter given by Eq. © is considered to be 77 « 5.26. 

Thus, by introducing the Forbes approximation for the SN barrier correction function v(f), the exponential factor 
in Eq. © can be expanded as follows: 


exp 


-u(f)b^ 2 /F 



exp 


—h(jA^ 2 /F 


(19) 
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FI G. 4 . Macroscopic current density, Jm (see Eq. as a function of the characteristic kernel current density, Jkc (see 

Eq.( 15 l), for 0.1 H ^ 0.9 and for an anode potential in the range 3.5kV < < lOfcF. The slopes (values of Jh - see Eq.(20l) 

for H = 0.1 and H = 0.9 were found to be 1.27 (dashed line) and 1.11 (solid line), respectively, with an error of 10 -3 in the 
linear regression fit (see Table II for other values of /3h)- The inset shows a snapshot of a portion of a LACFES with H = 0.1, 
in which the non-emitting locations (black) are distinguished from the emitting locations (white). The substantial increase in 
the effective area (see the text for more details) as the anode electric potential changes from <E>' 4 = 6 kV to 'F' 4 = 10 kV can be 
clearly observed. 


Fig. [4] shows the behavior of log 10 ( Jm) as a function of log 10 ( J^cO for 0.1 < H < 0.9 and an anode potential in the 
range 3.5 kV < $ A < 10 kV. These results suggest a scaling relation between Jm and Jkc as follows: 


Jm ~ [J k cf H {Ph > 1 ), ( 20 ) 

where the exponent fin depends on the H of the LACFES. According to Eq. [l4j this result can be understood as 
the result of a power-law dependence between a,w (equal to A m in our case) and Jkc, he., a M = Ci(Jfcc , ) C2 i where, 
ci > 0 and C 2 > 0 (ci = «« and C 2=0 if the emission area is constant) are constants in an appropriate range of the 
macroscopic electric fields. Indeed, if the surface of a LACFES exhibits fluctuations on small scales, then for a low 
electric potential of the anode (or, alternatively, a macroscopic electric field, Fm), only some points on the surface 
(those with larger FEF values) are emitters. If F M increases (decreases) by a factor A>1(0<A<1), J^c also 

increases (decreases) by a factor A 2-1 '/ 6 [exp (—60 3 / 2 /(tc-Fm))] Moreover, some new locations on the surface 

will become emitters in addition to those that were previously emitters. Thus, olm increases (decreases) by a factor 

/\( 2 -ri/ 6 )c 2 |g X p c 2 6^ 3 / 2 /(7cFm))] 1 1 • This is a very reasonable physical explanation because of the relatively 

slow decay which may be observed in the first and second regions of the local electric field distributions (for instance, 
corresponding to H = 0.1), which indicates that there are large regions on a LACFES where the differences in FEF 
values are relatively small. Thus, a small variation in the macroscopic electric field can result in the appearance of a 
non-linear relation between a« and Jkc- If the fluctuations of the FEFs on the LACFES are sufficiently large, then 
the emission locations with smaller FEF values may not be able to become field emitter, even at high experimental 
electric fields. The role of the exponent /3h (= 1 + C 2 ) is reflected in this behavior. 

It is evident that the linear dependence between Jm and Jkc (as a consequence of the lack of dependence between 
olm and Jkc, as shown in Eq. ( |14| ) ) is more approximate for 0.5 ^ H ^ 0.9, where the height fluctuations of 
the emitting surface are correlated (whereas they are random for H = 0.5). Outside this range, the presence of 
anticorrelations and, consequently, more small-scale irregularities of the LACFES contribute more strongly to the 
increase in the CFE area. The numerical calculations corroborate this assumption, and the corresponding values of 
/ 3h can be determined from the slopes of the curves shown in Fig. [4] The results indicate that for 0.1 ^ H < 0.5 
in particular, fin is significantly different from unity. The slopes for 0.1 ^ H ^ 0.9 are presented in Table II. This 
observation indicates that for rough LACFESs, small Hurst exponents play an important role in causing the rapid 
increase (relative to that observed in the case of smoother LACFESs) in the “effective area” across the surface of the 
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1/F m (nm/V) 


FIG. 5. Jm-Fm- type FN plot for LACFESs with 0.1 < H < 0.9. The corresponding slopes for H = 0.1 and H = 0.9 (dashed 
lines) are shown (see Table II for other values of Sf m {H) corresponding to 0.1 ^ H ^ 0.9). The work function of the LACFES 
is considered to be approximately constant and equal to (j> = 3.5 eV. The inset of this figure shows the derivatives of the 
Jm-Fm- type FN plots for LACFESs with Hurst exponents of H = 0.1 and H = 0.9. 


field emitter as the macroscopic electric field Fm (or, alternatively, <I>‘ 4 ) increases. This may be a signature of an 
interesting scale-invariant relation between Jm and Juc (as the result of a power-law dependence between om and 
Jfcc), which can motivate experimental tests with real morphologies that follow orthodox field emission. 

The inset of Fig. [4] shows a representation of a portion of a LACFES with H = 0.1, which illustrates the contrast 
between the non-emitting and emitting locations. To define these features, we apply a criterion such that a given 
location on the LACFES is an emitter if </&(</, F)/Jm > ICC 3 . The substantial increase in the effective area that 
occurs as the electric potential of the anode changes from = 6 kV to <F 4 = 10 kV is evident. This result is also a 
consequence of the two exponential regimes with small values of AS (see Table I) in the distribution p(y), which is 
more pronounced for 0.1 ^ H < 0.5, thereby increasing the probability of finding emitter locations with high FEF 
values. 

Now we address our primary purpose, namely, elucidating the effect of the Hurst exponent of an LACFES on the 
slope of the FN plot, which is typically used to extract relevant features of field emitter surfaces. Experimentally, 
the data obtained from field emission measurements can be described as orthodox when the characteristic field 
enhancement factor, yc, is independent of voltage. Indeed, in our model, this requirement is satisfied. Moreover, we 
assume that the emission is controlled solely by the tunneling barrier at the emitter/vacuum interface. Thus, our 
results can be tested experimentally under these conditions. Fig. [HJshows the behavior of In {(Jm/(Fm) 2 )} vs. 1 /Fm 
(or a Jm~Fm ~type FN plot) for a LACFES with 0.1 ^ H < 0.9. The FN plot appears to exhibit approximately linear 
behavior in the considered range of <E ,j4 , although in reality, the derivative is not constant, as observed in the inset of 
Fig. |5| for H = 0.1 and H = 0.9. An approximately constant value of the derivative is observed in the limit of low 
macroscopic electric fields, in which only the main peaks are emitters, followed by an increase in the corresponding 
absolute value in the limit of high values of Fm- 

The slope of the JM-^M-type FN plot (extracted from a linear least-squares fit to the numerical (or experimental) 
data), which depends on H, is denoted by Sf m (H). In fact, the higher is the characteristic FEF value, the lower in 
magnitude the slope of the FN plot will be. This trend is evident in Fig. [5] and reflects the effect of the different 
Hurst exponents of the LACFESs (for the same global roughness {W ~ 75 nm) - S(f m (0.5) > S Fm (0.3) > S Fm (0.2) > 
S(f m ( 0.1)), thereby corroborating the inequality in the corresponding yc values shown in Table I. The values of 
Sf m (H) are presented in Table II. 

First, we use the elementary slope characterization parameter (SCP), which is denoted by y<£ (see Ref.[4]), to 
simulate the common practice of experimentalists in using the triangular-barrier (or elementary) FN equation (Eqs. 
([8]) and Q for A cfe = Ef = 1). In this case, the SCP is given by 


lc = _ b4?/ 2 /S FM {H ). 


(21) 
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H 

S FM (H)(V/nm) 

T 

7 c 

Ph 

Ph ,cr 

7 c 

0.1 

-4.60 

± 

0.05 

9.72 

± 

0.05 

1.271 

± 

0.001 

11.79 ±0.05 

0.2 

-4.88 

± 

0.03 

9.16 

± 

0.03 

1.185 

± 

0.001 

10.32 ±0.03 

0.3 

-5.32 

± 

0.02 

8.40 

± 

0.01 

1.174 

± 

0.001 

9.39 ±0.01 

0.4 

-5.45 

± 

0.02 

8.20 

± 

0.01 

1.132 

± 

0.001 

8.87 ±0.01 

0.5 

-5.60 

± 

0.01 

7.98 

± 

0.01 

1.127 

± 

0.001 

8.59 ±0.01 

0.6 

-5.87 

± 

0.02 

7.61 

± 

0.02 

1.145 

± 

0.001 

8.33 ±0.02 

0.7 

-5.92 

± 

0.02 

7.55 

± 

0.02 

1.124 

± 

0.001 

8.10 ±0.02 

0.9 

-6.06 

± 

0.01 

7.37 

± 

0.01 

1.113 

± 

0.001 

7.82 ±0.01 


TABLE II. Results for 0.1 ^ H ^ 0.9. The slopes of the J m-Fm -type FN plots shown in Fig. [5| the characteristic FE F v alues 
7 cr, the /3h values extracted from the results presented in FigUland 7c and 7c H <T were calculated using Eqs. |2l| and 

( |23| ) while considering the elementary FN equation and relevant corrections, including the SN barrier and the morphology of 
the LACFES, respectively. 


From the slopes presented in Fig. [ 5 ] (see Table II) and using Eq. (211, we calculate 7 ^ = 7.37 (error of 5% compared 
with the 7 c value reported in Table I) for H = 0.9 and 7 ^ = 9.72 (error of ss 20% compared with the 7<7 value 
reported in Table I) for H = 0.1. This latter error, which is related to the fact that there is a greater probability of 
finding field emitter locations on the LACFES with larger FEF values than is naively expected (as is also the case for 
H values in the interval 0.1 < H < 0.4), represents an underestimation of the real characteristic FEF value that is 
of practical significance. Table II summarizes the values of 7 ^ calculated for 0.1 ^ H ^ 0.9. To further explain this 
result, let us return to the form of the scaling between Jm and Jkc■ From Eqs. § and p0||, it follows that 


ln{(J M /(E M ) 2 )} = r 


d (in [J m (H)/Fm]) 2 (/ 3 if-i) 


d ( 1 /E m ) 

from which it is possible to define the corrected FEF: 


F 

5 r M 


f /3 H v Fc b(t> 3/2 \ 1 

l 7 tr" r*M' 


( 22 ) 


ft H ,<T 

7c 


/3 H vb(j) 3/2 \ 
Sf m (H) ) 


(23) 


In Eq. (23), a is a generalized slope correction factor related to the SN barrier and /3 h includes the effect of the 


geometry of the LACFES on the estimation of the characteristic FEF. Formally, the first term in Eq. (22) does not 
define the intercept in the FN plot. This occurs only if a triangular barrier, which corresponds to an unrealistic 
physical situation, is considered and if the area of emission of the rough LACFES is independent of voltage. This is 
not the case in our system, where the effects that can produce non-linear features in the .//vf-AW-type FN plot are 
the dependence between the area of emission of the rough LACFES and the anode voltage as well as the effect of the 
SN barrier. According to this argument, the macroscopic current density is a function of H such that Jm = 

The process of deriving emission-area estimates from FN plots has previously been investigated by Forbes et al. [53] 
using free-electron theory and considering three different tunneling-barrier models. However, the current work can 
motivate the investigation of an additional effect, namely, that of the irregular morphology of the LACFES. In this 
way, it is possible to define, based on the function T in Eq. (22), an effective Area Efficiency of Emission, a e ff, 

that clearly depends on both the macroscopic electric field, Fm, and the derivative a ( ln ^A^/ F M]) . This quantity, 
cteff, can be extracted from the term F ^ ^ , F^f H ~ 1 ^^ in a given range of values of the variable 1 /Fm, 

again because most emitters melt for values of 1 /Fm less than some reference value [34j (and this effect may be more 
pronounced in the case of a rough LACFES). 

Obviously, because /3h is close to unity (F' 2 j' 3ff_1 ' ) — » 1), a e ff ~ Interestingly, this feature also appears, 

more approximately, in the case of 0.5 ^ H ^ 0.9, in which the height fluctuations on the LACFES are statistically 
correlated (or random), a typically takes a mid-operating-range value of approximately 0.955 0134]. If we apply 
this assumption (which is certainly a good approxi mati on for our system), consider the values of the FN plot slopes 
obtained from Fig. [ 5 ] (see Table II) and use Eq. (23) to estimate the corrected characteristic FEF value, , 
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then we obtain the values presented in Table II. Impressive agreement is clearly evident between r )£ H,a and the 
values presented in Table I that were extracted from the distribution p( 7 ), which is not accessible to field emission 
experimentalists. Thus, the correction (3h must be applied (primarily for LACFESs with large amounts of individual 
field emitter tips) to more precisely extract the characteristic FEF of LACFESs from FN plots. The effect discussed 
above can, of course, alternatively be described as being caused by the variation in the effective emission area with 
voltage (or, equivalently, with the macroscopic field), and our conclusion is that this variation must be taken into 
account. 

Alternatively, when empirical CFE (i-V) characteristics are used (where “V” is the measured voltage) , our 
results suggest replacing k = 2 — p((j))/6 (the result from Eq. ( fl7| ) , considering the assumptions that lead to Eqs. 
(14) and (15) as well as the expansion with k(H) = (2 — r](^)/6)pH for a rough LACFES (for constant /3h), if 
the range of the macroscopic electric field yields an “adequately linear” FN plot. This suggests that an investigation 
of the real irregular morphology of the LACFES, through the precise measurement of the Hurst exponent (and the 
global roughness), is an important method of gaining a more complete understanding of the relevant measurements to 
be extracted from the FN plot to evaluate the electron emission capability of a material with an irregular morphology. 
Moreover, a systematic and careful investigation of the scaling between Jm and JkCi which are both experimentally 
accessible measurements, is of fundamental importance for evaluating the effects of the rough geometry of the LACFES 
on the corresponding CFE properties. Situations in which olm scales as a more complicated function of Jkc cannot 
be disregarded, and an investigation of these features in potential CFE materials is certainly desirable. Finally, our 
model can also be applied to investigate the inclusion of work function distributions over rough LACFESs, which has 
not been considered here. The investigation of this effect will be a subject of future work. 

In summary, we investigated the role of the Hurst exponents of rough LACFESs on the corresponding electric field 
distributions and the field emission quantities of technological interest, particularly the characteristic FEF value, 7 c- 
Our results, which were obtained by considering the orthodox CFE and LACFESs with a global roughness on the order 
of a few tens of nanometers, demonstrated that correlated (or random) LACFES morphologies (0.5 ^ H ^ 0.9) exhibit 
a more weak dependence between the area of emission and the macroscopic electric field, whereas for 0.1 ^ H < 0.5, 
a stronger dependence is evident that contributes to the formation of non-linear features in the FN plot. This is a 
consequence of small-scale fluctuations in the morphology of the field emitters, even when all LACFESs have the same 
global roughness. For any H 1 the local electric field distribution over the LACFES exhibits two exponential regimes 
and a power-law scaling between the macroscopic current density ( Jm) and the characteristic kernel current density 
( Jkc ) f° r the typical experimental range of macroscopic electric fields. Our results indicate a scaling of the form 
Jm ~ [Jkcf H , where the exponent fin > 1 depends on H. Moreover, for orthodox field emission from a LACFES, 
this scale must be considered when calculating the slope characterization parameter (in addition to the generalized 
slope correction from the SN barrier) that is used by experimentalists to extract the characteristic FEF. Failing to 
include this parameter in the SCP (using FN plot data) may result, particularly for LACFESs with 0.1 ^ H < 0.5, 
in a significant underestimation of the characteristic FEF, which is a very useful measure of the emitter sharpness. 
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Supplementary Information I 
Numerical method for electric field calculations 

Assuming that the region between the emitter surface and the anode is vacuum, this implies the numerical so¬ 
lution of Laplace’s equation, 


V 2 $ = 0, (24) 

in a discretized space, imposing suitable boundary conditions. In our case, they correspond to Dirichlet conditions 

both at the cathode (<1> = 0) and the anode ($ = < 1> A ), and periodic lateral conditions, so that for a domain formed 
by L x x L y x L z points (in the directions a ~,y,z, respectively), we have 


— *&i+L x ,j,k — *&i,j+L y) k 


/j x ,j + Ly,k- 


(25) 


Equation (24) is then solved iteratively using a second order finite difference scheme, in which the electric potential 
is given by 
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That is, at each iteration the new potential is just the average of its magnitude on the surrounding points at 

the previous iteration, (to). These values are then refined iteratively starting from the first step, = 0, until a 

satisfactory convergence criterion 
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is met. In this work, the error is set to be £ = 10 . In Eq.(26), ui is the length of lattice parameter in the cubic 


lattice. In this work, we refine the grid in such that u> = 1/50 u, where u is the basic unit distance. However, no 
significatively deviations has been detected using u> = u. Once convergence has been achieved, the intensity of the 
electric field, 


F(r) = -V$(r), (29) 

at a any point (not necessarily on the grid), r p = (xp,yp, zp), can be evaluated by linear interpolation from the 

values at the eight vertices of the circumscribing grid cube. If such a cube “starts” at the grid location ( i,j,k ), the 
corresponding cartesian components, F x,y,z { r), are given by 


F»(r) =N 


uj rpH 

^i+aij+aj ,k+a k 


(30) 
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FIG. 6. Characteristic FEF, evaluated numerically by the solution of the Laplace equation, as a function of the anode electric 
potential, considering a semi-sphere and a cone protuberances superimposed to a planar conducting surface. 
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and 


1"(xp, ai ,y P , aj ,ZP,a k ) = \J(Xp - x i+ai ) 1 2 + (yp - Uj+a-j) 2 + ( Z P - Z k+ak ) 2 


(32) 


To test the accuracy of our method, we consider calculate the characteristic FEF, 7 c, solving the electric field 
in three dimensions around a conducting sphere (and conical) structure superimposed in the center of a planar 
conducting substrate. The electric potential at the surface of the conducting sphere was fixed to $ = 0, while at the 
anode a electric potential <F A . The results for 7 c of the two structures are presented in Fig. [ 6 j for Axi,Axj -C T, 
as a function of the electric potential $ A . It’s clear the very good agreement between our numerical results and the 
analytical solution which predicts 7 c = 3 for the semi-sphere on a plane |T. In the case of conical structure, our 
results are in accordance with those presented on reference [ 2 ]. 
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Considerations regarding three dimensional rough LACFESs 
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In this Supplementary Information section, we wish to demonstrate that the same relevant results obtained in 
the one-dimensional case (using orthodox cold field electron emission theory) are also obtained for a genuine rough 
two-dimensional LACFES. For this purpose, we use fractional Brownian motion (fBm) algorithms jT] to generate 
random self-affine objects with specific Hurst exponents H. FBm algorithms can be generalized to higher dimensions 
using a multidimensional process such that for three dimensions, an irregular LACFES of height h(x 7 y) results in 
interfaces that are statistically equivalent in all directions. We consider the self-affine sets that are constructed using 
the well-known midpoint displacement algorithm for H= 0.1 [2J. 

In Fig. [TJ we illustrate the rough LACFESs produced using the described fBm algorithm for a Hurst exponent 
of H=0.1 and with global roughnesses W « 20 ran, 37.5 nm and 0.5/im. In Fig. |8j several equipotential surfaces 
are shown for the case in which the electric potentials of the LACFES (with H = 0.1 and W ~ 37.5 nm) and the 
far-away conducting anode are < F‘ ? = 0 V and <I>' 4 = 500 V, respectively, corresponding to Dirichlet conditions. We 
apply this approximation, following the same procedure used in the two-dimensional case, with the intent of rounding 
the sharpest projections. In this way, we replace the real fBm surface with an equipotential surface, <f> £ = 1 V (very 
similar to the previous one), that represents an approximation of the field emitter surface. Again, following the same 
procedure using the WM function, this methodology ensures that on small scales (though larger than the atomic 
scale), the local roughness is negligible, whereas on large scales (though much less than the lateral size of the system, 
L), the local roughness scales with the same exponent H. 

In Fig. |9j we present the local FEF distributions, for LACFESs shown in Fig. [7] Interestingly, exponential 
behaviors (two exponentially decaying regions as well as a characteristic FEF factor of the same order, for W ~ 37.5 
nm and W ~ 0.5/mr, as that observed in thin-film emitters with irregular surfaces - see Ref. [30] of the manuscript) are 
observed in all cases, the same behavior as that discussed in the case of WM surfaces, when the one-dimensional field 
variation across the emitter surface is considered. Thus, the one-dimensional case is an interesting tool that captures 
the main results experimentally obtained for rough LACFESs, as noted previously. 

Note that for H = 0.1 and W ~ 37.5 nm [« 0.5/nn] (in which case the characteristic FEF has been found to be 
7 c « 97.24 [« 519.74]), two exponentially decaying regimes are evident, labeled as (I) and (II), which correspond to 
p T ( 7 ) ~ exp(—S^j) and p JI ( 7 ) ~ exp(-fijfy), respectively, where > 5 ! H . For W ~ 37.5 [« 0.5/xm], Sq 1 xlog 10 (e) = 
(0.0324±0.0001) [<5 q j x log 10 (e) = (0.00273±0.00005)] is identified for interval (I), and S& x log 10 (e) = (0.046A0.001) 
[6 q\ x log 10 (e) = (0.0090 ± 0.0001)] is found for interval (II). By contrast, for W « 20 nm, only one exponentially 
decaying regime is observed, with 6 ^ A x log 10 (e) = (0.247 ± 0.008). This finding indicates that the p(y) distribution 
is also sensitive to changes in the global roughness of the LACFES. Differences between the theoretically obtained 
values of the characteristic FEF and those obtained experimentally can also be caused by differences in the global 
roughness. Table III summarizes the values of the parameters extracted from the p(y) distributions presented in Fig. 
[9] and the 7 c values for all values of W explored. 

Finally, we show that the fiw (note that the notation is now used with W for global roughness) correction, which 
includes the effect of the geometry of the three-dimensional LACFES on the estimation of the characteristic FEF, 
must be considered to ensure more precise estimation of 7 c (which was found to be 7 c ~ 97.24 [ss 519.74] from the 
p( 7 ) distribution depicted in Fig. [ 9 ] for W ~ 37.5 nm [ss 0.5/Ltm]). We restrict this discussion to H = 0.1 and for 
roughnesses W ~ 37.5 nm and W ~ 0.5^m, i.e., to the LACFESs represented in Figs. [7] (b) and (c), respectively. 
Fig. [Hjja), shows the behavior of the corresponding Jm-Fm~ type FN plots. The FN plots seems, again, to exhibit 
approximately linear behavior for the considered range of macroscopic electric field Fm, although, in reality, the 
derivative is not constant, as observed in the corresponding inset. 

Fig. 10 ]b) shows the behavior of log 10 (Jm) as a function of log 10 (Jfcc), considering a macroscopic electric field 
in the range 0.07 V/nm < Fm < 0.1 V/nm and 12.5 V/pm < Fm < 2QV/pm. for W ~ 37.5 nm and W « 0.5^m, 
respectively. These results suggest, again, a scaling relation between Jm and Jkc with /3w = (1-255 ± 0.001) for 
W ~ 37.5 nm [and ftw = (1-222 ± 0.002) for W ~ 0.5/rm]. From these results, for W ~ 37.5 nm [« 0.5/um] we 
calculate an elementary slope characterization parameter of 7 ^ « 77.7 [« 434] (error of 21% [16.4 %] with respect to 
the 7 c value found from the p( 7 ) distribution). If we include the corrections a (which is a generalized slope correction 
factor related to the SN barrier) and /3w (which includes the effect of the global roughness of the LACFES on the 
estimation of the characteristic FEF), then the corrected slope characterization parameter is found to be 7 ^A’ <T ~ 93 
[« 505.8], yielding an error of 4% [2.6%] with respect to the 7 c value determined from the p( 7 ) distribution (see Table 


IV for values of the slopes of the Jm-Fm -type FN plots shown in Fig. 
Pw values extracted from the results presented in Fig 


10 


10 

(b) and 7 ^). 


(a), the characteristic FEF values 7 ^, the 


In conclusion, the scale that results in the exponent ftw must be considered (including in the case of a three- 
dimensional setup and in addition to the generalized slope correction from the SN barrier) in the calculation of the 
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w 

&h x log 10 (e) 

$h x l°g 10 (e) 

7 c 

~ 20 nm 

0.247 ±0.008 

- 

w 38.35 

« 37.5 nm 

0.0324 ±0.0001 

0.046 ±0.001 

« 97.24 

« 0.5pm 

0.00273 ± 0.00005 

0.0090 ± 0.0001 

« 519.74 


TABLE III. Results of the extraction of parameters from the p( 7 ) distribution shown on Fig. [ 9 ] for H = 0.1 and global 
roughnesses W « 20 nm, 37.5 nm and 0.5 pm. For W « 20 nm a dominant exponential decay can be observed in the p( 7 ) 
distribution, characterized by 5^. 


W 

S Fu (W)(V/nm) 

T 

7 c 

/3w 

7 

« 37.5 nm 

-0.5759 ±0.0009 

77.77 ± 0.01 

1.255 ±0.001 

93.2 ±0.3 

~ 0.5pm 

-0.1031 ±0.002 

434.2 ± 0.1 

1.224 ±0.001 

505.88 ± 0.01 


TABLE IV. Results for H = 0.1 and W 

Figliol(a) 


20 nm, 37.5 nm and 0.5pm. The slopes of the type FN plots shown on 


characteristic FEF yj, the /3w values extracted from the results presented on Fig. 


10 


b) and 7 ^ W,CT . 7 c an( i 7 c w ’ 

were calculated using Eqs. (21) and (23) while considering the elementary FN equation and relevant corrections, including the 
SN barrier and the morphology of the LACFES, respectively. 


slope characterization parameter that is used by experimentalists to extract a more precise characteristic FEF, such 
as that found using the one-dimensional WM function. 
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FIG. 7. Portion of rough LACFES generated using the fBm algorithm for a Hurst exponent of H=0.1 and global roughnesses: 
(a) W « 20 nm; (b) W ss 37.5 nm and (c) W « 0.5/rm. 
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FIG. 8. (top) Snapshot showing several equipotential surfaces calculated numerically from the solution of the Laplace equation, 
using the LACFES represented in Fig. [t] (b) with appropriate Dirichlet conditions (<I> S = 0 V and <F A = 500 V - see the text 
and Supplementary Information I for more details). The color bars indicate the electric potential values (in Volts), (bottom) 
The equipotential surface that is used as an approximation to the field emitter surface corresponding to the real fBm surface, 
defined by & E = 1 V. The color bars indicate the height of LACFES (in nanometers). 
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FIG. 9. Local FEF distributions, p( 7 ), for LACFESs constructed using the fBm algorithm for H = 0.1 with global roughnesses 
of W fa 20 nm (red squares), W fa 37.5 nm (blue triangles) and W fa 0.5/rm (green hexagons). Exponentially decaying behavior 
can be clearly observed in all cases. In particular, for W fa 37.5 nm and W fa 0.5/rm, the same behavior previously observed, 
that of two exponentially decaying regimes, is apparent. The slope of the dashed (red) line is (—0.247 ± 0.008). The slopes 
of the solid (blue) lines are (—0.0324 ± 0.0001) (region (I)) and (—0.046 ± 0.001) (region (II)). The slopes of the dot-dashed 
(green) lines are (—0.00273 ± 0.00005) (region (I)) and (0.0090 ± 0.0001) (region (II)) (see Table III). The differences observed 
in the p(y) distributions are essentially related to the difference in global roughness between the LACFESs. 
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FIG. 10. (a) Jm-Fm- type FN plot for LACFESs with H = 0.1 and W « 37.5 run (blue triangles) and VF « 0.5nnr (green 
hexagons). The corresponding slopes (dashed lines) are shown (see Table IV). The work function of the LACFES is considered 
to be approximately constant and equal to <f> = 3.5 eV. The inset of this figure shows the corresponding derivatives of the 
Jm-Fm- type FN plots for LACFESs. (b) Macroscopic current density, Jm as a function of the characteristic kernel current 
density, Jkc with H = 0.1 and W ta 37.5 nm (blue triangles) and W ~ 0.5/rm (green hexagons). For W « 37.5 nm and 
W « 0.5/rm the macroscopic electric fields are in the range 0.07 V/nm < Fm < 0.1 V/nm and 12.5V/[im < Fm < 20V//rm, 
respectively. The slopes (values of /3w ) ar e also indicated (see Table IV). 












































